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Abstract: We present a method for the computation of hepta-cuts of two loop scattering 
amphtudes. Four dimensional unitarity cuts are used to factorise the integrand onto the 
product of six tree-level amplitudes evaluated at complex momentum values. Using Gram 
matrix constraints we derive a general parameterisation of the integrand which can be 
computed using polynomial fitting techniques. The resulting expression is further reduced 
to master integrals using conventional integration by parts methods. We consider both 
planar and non-planar topologies for 2—7-2 scattering processes and apply the method 
to compute hepta-cut contributions to gluon-gluon scattering in Yang-Mills theory with 
adjoint fermions and scalars. 
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1 Introduction 

Precision cross section predictions for hadron colliders are an essential tool in the search for 
new physics. While one-loop amplitudes give access to quantitative Next-to-Leading Order 
(NLO) background estimates, a reliable analysis of the theoretical uncertainty requires 
Next-to-Next-to-Leading Order (NNLO) corrections. 

The use of Feynman diagram techniques for the evaluation of scattering amplitudes has 
always presented a major challenge owing to the rapid growth in complexity with increasing 
loop order and external legs. In recent years the development of on-shell methods [1-5] has 
played a major role in removing this traditional bottleneck for both tree-level and one-loop 
amplitudes.^ Generalised unitarity [3, 9-18] and integrand reduction techniques (OPP) 
[19] have been developed into fully automated numerical algorithms able to compute high 
multiplicity one-loop amplitudes [20-29]. 

Unitarity methods for multi-loop amplitudes have proven to be extremely powerful 
tools in super-symmetric gauge theories. Maximal cutting techniques are an efficient 
method for reducing these complicated amplitudes to the evaluation of a limited num- 
ber of master integrals. These techniques have been developed in the course of gluon-gluon 
scattering amplitudes in A/" = 4 Super- Yang-Mills (SYM) enabling computations up to five 
loops [30-33].^ At two loops computations with up to six external legs have been achieved 
[34-36]. Octa-cuts [37] and the related leading singularity method [38, 39] are also valid 
approaches in = 4 SYM enabling the computation of loop amplitudes directly from 
tree- level input. 

In non super-symmetric theories, like QCD, the basis of integrals is far more compli- 
cated yet the current state-of-the art techniques have been able to compute 2 — )• 2 processes 
in massless QCD [33, 40-45]. The motivation for the present study is to use some of the 
technology successful in the super-symmetric cases to simplify the computation of these 
amplitudes. 

The aim is to construct full two-loop amplitudes from products of tree-level amplitudes 
following the successful approach taken at one- loop. Following a top down approach one 
begins with the leading singularities, then systematically reduces the number of cuts to 
study more of the full amplitude. At each step one subtracts the singularity structure 
previously constructed in order to obtain a polynomial system. At one-loop this procedure 
relies on the knowledge of a basis of integral functions. Though such a basis is not known 
at two loops, the reduction of arbitrary loop integrals can be understood using integration 
by parts (IBP) relations [46, 47]. Using a restricted set of IBPs constructed using Grobner 
bases, Gluza, Kajda and Kosower were able to construct a unitarity compatible integral 

^See [6-8] and references therein for recent reviews on these topics. 

^We note that only the coefficients of the master integrals are known at five loops, not the integrals 
themselves. 
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basis for planar topologies [48]. Schabinger recently showed similar sets of IBPs could be 
obtained without the use of Grobner bases [49] . To date both a maximal unitarity approach 
[50] and an integrand reduction program similar to OPP have been proposed [51] which 
explore the use of fitting such a basis from tree-level input. 

The approach we follow here will allow us to construct a general integrand parame- 
terisation from analysis of Gram matrices. This system can be matched to an expansion 
of the products of tree-level amplitudes evaluated at a complete set of on-shell solutions 
to the loop momenta. This leads to a linear system of equations that can be inverted to 
derive a master formula for the reduction of the integrand. The part of the integrand that 
remains after integration is compatible with further reduction to master integrals by using 
conventional IBP identities and Lorentz Invariance identities [52] by means of the Laporta 
algorithm [53], or the related approaches [54-56]. A number of public tools [57-60] are 
available to perform this step of the computation. 

We address the first in a long list of ingredients required for a general decomposition 
of two loop amplitudes, the maximum singularities in 2 — )• 2 processes. These are all 
contributions with seven propagators that can be extracted from hepta-cuts in four dimen- 
sions. The procedure reduces the computation of the amplitude to a polynomial fitting 
procedure over a product of tree-level amplitudes and is amenable to both analytic and 
numerical techniques. We compare with the super-symmetric results obtained using the 
recent approach of Kosower and Larsen [50] . 

Our paper is organised as follows. We begin by re-deriving some results of the gen- 
eralised unitarity algorithm at one-loop. We focus on some of the key issues that we will 
apply to the two-loop case. We then turn out attention to the three independent seven 
propagator topologies for 2—7-2 processes, the planar double box and penta-box config- 
urations as well as the non-planar crossed box configuration. We develop a method for 
determination of an integrand parameterisation using constraints from 5x5 Gram matri- 
ces. We then use this information to construct an invertible linear system from the full set 
of on-shell solutions which maps the this parameterisation to products of tree-level ampli- 
tudes. The resulting integrand can be further reduced to master integrals by application 
of well known integration by parts identities. We demonstrate the technique by applying it 
to the four-gluon scattering in Yang-Mills theory. The expressions can be related to those 
in super-symmetric Yang-Mills and we comment on some the simplifications that occur in 
those cases. Finally we present our conclusions and some outlook for future studies. 

2 Review of Generalised Unitarity at One-Loop 

In this section we will re-derive the well known integrand parameterisation used in numeri- 
cal one-loop generalised unitarity algorithms [9] and the closely related integrand reduction 
of Ossola, Papadopoulos and Pittau (OPP) [19]. For more detailed reviews of the subject 
we refer the reader to refs. [6, 7]. 

We represent a general ordered one-loop amplitude as a product of rational coefficients 
multiplying scalar integral functions with four or fewer propagators. For the present exer- 
cise we will restrict ourselves to cases where all numerators are in four dimensions and the 
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pentagon contributions are collected into a remaining rational contribution, 
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+ tadpoles, wave-function bubbles and rational terms. 



(2.1) 



In the above we have defined the inverse propagators Di^ = {k—pi^^i^-i)"^ and the dimension 
d = 4 — 2e. We define pij = X^'^^j Pk as the sum of external momenta such that pi-^/i_^^\ = 
and have taken the restriction that all propagators are massless. 

A general one-loop amplitude can be computed by repeated evaluations of a process 
specific numerator once general forms for the integrands, Ac,x(fc)) have been constructed. 
The numerator used to fit the cut integrands could be generated from a Feynman diagram 
representation but in the following we will take a top down approach and factorise each 
cut into products of tree- level amplitudes. 

We will go through this known procedure in some detail since our generalisation to 
two loops will follow it closely. 



2.1 Quadruple Cuts 



Pi 



P3 



P2 



Figure 1. Conventions for the momentum flow in the one- loop box. 

Quadruple cuts of the one-loop amplitudes were first considered in the work of Britto, 
Cachazo and Feng [3]. Each integrand, A4^j^j2j3«4 5 depends upon three independent external 
momenta, say {Pi =pi^^i^-i,P2 = pi^^i^^i, P'i = Pi2,ia~i] as indicated in fig. 1. 

In order to span the full four dimensional space of the integrand we are able to define 
a vector w, satisfying w • Pfc = and > 0. Such a direction can be called spurious since, 
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{27T)d Di,Di,D,,Di, 



0. 



(2.2) 
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A simple representation for this complex vector is given by the totally anti-symmetric 
tensor w'^ oc e'^^^^. Equivalently, using the basis of massless vectors (i^Jj-ftTg) constructed 
from Pi and P2 as described in Appendix A, we can write: 

^^{P^) = ^ [{K\\P^\K\\{K\mK\\ - {K\\P^\K\\{K\mK\^ . (2.3) 



2712 

This gives a set of scalar products with which we can write down a completely general form 
of the cut integrand, 

{{k ■ k), (k ■ Pi), (fc • P2), {k ■ P3), {k ■ io)}. (2.4) 

We are able to re-write some of the dot products in terms of inverse propagators and 
constant factors, 

{k-k) = D,„ (2.5) 

2{k ■ Pi) = A4 - D^^ - Pi, (2.6) 

2{k-P2) = D,,-D,,+Pi, (2.7) 

2{k ■ P3) = Di, - As + 2P2 • P3 + P|, (2.8) 

where all Di^ vanish when the on-shell conditions are applied, 

{e = 0,{k- P^f = 0, {k - P2 - Psf = 0,ik + Pif = 0}. (2.9) 

This leaves us with one irreducible scalar product (ISP), 

a 

Renormalizability tells us that q < 4, yet we must find another relation before we are in a 
position to apply the cuts. Such information can be simply extracted from Gram matrices 
[9]. For 2n vectors {/i, ...,/„; vi, ... , Vn}, the n x n Gram matrix G is defined as 

G^cf ••• ^" ) , G^,=k-vj. (2.11) 

yvi ... Vn J 

In particular, for the case where {h, . . . , In} is identical to {vi, . . . , Vn}, we define 

G{h,...,ln)^G(\'---'\A. (2.12) 

y '1, • • ■ , 'n J 

The determinant det G is linear and anti-symmetric in the vectors in each row, 

detCT^ +detGri ■■■ in\ (2.13) 

yvi ... Vn J yvi ... Vn J 

-detG(^' ■■■ ^".^ (2.14) 

\ViV2...Vn / 
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( h + I'l . 


■ ■ 




\ Vl . 
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Therefore det G vanishes if {/i, . . . , /„} or {vi, . . . , Vn} are hnearly dependent. 

A Gram matrix can be used to calculate the inner products of two vectors, if their 
projection on a given basis is known. Let {ei,...,ed} span a d-dimensional vector space 
and let / and v be two vectors in that space. Once we expand v = ViCi and I = liCi and 
define Gd=G(ei, ...,6^) we can write, 



fv 


ei\ 
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(2.15) 



and so. 



(/•7;) = (/i,, 



, ld)G, 




{l-ei,...,l-ed)G,^ 




(2.16) 



If the basis is orthonormal, the Gram matrix becomes the identity and the relation 
eq.(2.16) is trivial. 

Explicitly, let d = 4 and {Pi, P2, P^jUj} be the basis of 4-dimensional loop momenta. 
When k is 4-dimensional momenta, by eq.(2.16), we obtain a quadratic relation of the 
Lorentz invariants, 



{k ■ Pi,k ■ P2,k ■ P^,k ■ uo)Gl^{k -Puk- P2, k-P3,k- Lof. 



(2.17) 



Alternatively, when /c is a 4-dimensional loop momenta, by the linear dependence property, 



detG 



Pi P2P3Ujk 

Pi P2P3^h 



0. 



(2.18) 



It is easy to see that (2.18), is equivalent to (2.17). We can get a non-trivial relation by 
computing the on-shell constraint fc^ = 0, 

{k-ujf = -u}'^{k ■ Pi,k ■ P2,k ■ P3)G^^{k ■ Pi,k ■ P2,k ■ Psf = const. (2.19) 

where G3 = G(Pi,P2)^3) is a constant matrix and {k ■ Pi), i = 1,2,3 are all constant at 
the quadruple cut. This relation tells us the a = 0, 1 in eq. (2.10). 

This is the maximum number of constraints available for this topology so we can turn 
our attention to the on-shell solutions for the loop momentum which will allow us to fit 
the coefficients cq and ci, the only coefficients left in (2.10) by the constraints. 

Following the well known parameterisation from the literature in terms of two-component 
Weyl spinors we find two complex solutions for k satisfying the constraints of eq. (2.9). On 
each solutions, of which explicit forms are given in Appendix C, the amplitude factorises 
onto a product of tree amplitudes, Ti^jaisM! 



(A;W) = CO + ci(fc« • oj) = Ti,i,i,i,ik^'^) = ds. 



(2.20) 
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where, 



Afe=± k=l 

(2.21) 

with n + 1 = 1. Since the four on-shell constraints, eq. (2.9), freeze the loop momentum, dg 
are just complex numbers. Prom the explicit solutions (see Appendix C) we find (A;«-a;) = 
\/V4 and (fc(^^ • oj) = —\fVi which, after feeding into eq. (2.20), leads quickly to the final 
result: 

CO = ^ (di + d2) , (2.22) 
ci = ^ (di - da) . (2.23) 

V M 



The important feature of this analysis was the use of the Gram matrix to constrain the form 
of the integrand which was then mapped to the products of tree amplitudes via the on- 
shell solutions to the loop momentum. For this simple case the number of loop momentum 
solutions and the number of independent coefficients in the integrand were the same. As 
we will see in the next example this feature is not always true. 

2.2 Triple Cuts 

To re-derive the formula for the triple cut integrands, Aa^j^jjis) '^^ follow exactly the same 
procedure as above. This time our space is spanned by two external momenta and two 
trivial space vectors v = {Pi = pj3^j^_i,P2 = Pii,i2-i^^i:^2} where in the {K'1,K2) basis 
we can write, 

^'i=l {{K\\r\K'2] + {Kh'\K\]) , (2.24) 



{K\\j'^\kI]-{K',\j'^\K\]). (2.25) 



Removing the trivial scalar products that can be re-written in terms of propagators we 
have the following form for the integrand, 

A3,,,i,i,{k) = Y,Cai3{k-u;ir{k-u;2f. (2.26) 

Renormalizability implies that a + /3 < 3 and therefore there are ten Cij coefficients. The 
Gram matrix identity, 

\ Fi F2 uji UJ2 k j 

leads us to the relation, 

{k ■ ujif + {k ■ uj2f = const. (2.28) 
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which reduces the number of independent Cq,^ coefficients to seven. In principle we could 
chose any seven of the ten but in order to ensure terms in the integrand proportional to 
{k ■ Uk) integrate to zero, it is convenient to choose: 

C = (coo, Coi, Cio, C20;02, Cll, Ci2, C21), (2.29) 
so that integrand is written, 

^3,ili2iaif^) = ^00 + Clo{k ■ UJi) + Coi(/c • W2) + Cii{k ■ UJi){k ■ UJ2) 
+ Ci2{k ■ U}i){k ■ UJ2f + C2l(/c • Ulif{k ■ IJJ2) 

+ C20;02[{k-UOif -{k-UJ2?) . (2.30) 

We then turn to the on-shell constraints, 

{k" = 0,{k- Pi? = 0,{k-Pi- P2? = 0}, (2.31) 

of which there are two possible solutions: 

= aiKj'^ + a2K'^' + \{K\mK'2] + (2.32) 
= a,KY + a2KY + {{K^mKW + ^{K\mKll (2.33) 

where, 

By feeding this into eq.(2.30) we define the coefficients d^^x which can be extracted from 
the subtracted product of three tree-level amplitudes, 



x=—3 



Equating coefficients of t on both sides of this equation gives us a 14 x 7 matrix, M, which 
relates the dg^x coefficients to the coefficients, 

d = M • c (2.36) 

where d = ((ii _3, di _2, di _i, di^o, £^1,1,^1,2, c?i,3, ^2 -3, ^2 -2, ^2,-1,^2,0, ^^2,1,^2,2, ^2,3)- The 
final step is to invert M, the fact that it is invertible means that we have a unitarity cut 
compatible parameterisation of the integrand. The inverse falls into two regions, firstly 
when all P| 7^ (corresponding to 0102 7^ 0) the null space of M contains all of the k^'^^ 
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solution and c coefficients are, 



Coo — difi 

ClO = di i - 0102^1 3 —di I 

aia2 a 

coi = idi^i + iaia2di 3 H di i + 

aia2 

2i , 

cii = 2zdi,2 + -ir^di-2 

C20;02 = (^1,2 2~2'^l,-2 

«I«2 

4i 

C21 = 4zdi,3 + ^rT*^! -3 

4 

Ci2 = -4(ii 3 + ^^di _3. 



2 2'^l'-3 


(2.37) 


(2.38) 


'l'^2 




2 2'^l:-3 


(2.39) 


1 Z 






(2.40) 




(2.41) 




(2.42) 




(2.43) 



Coo 


= di,o 


cio 


= - ^2,1 


coi 


= i ((ii,i + (i2,l) 


Cll 


= 2i (di,2 + (i2,2) 


C20;02 


= (ii,2 — ^2,2 


C21 


= 4z ((ii,3 + (^2,3) 


Cl2 


= -C?l,3 + C?2,3 



In case any = we will find 0102 = and all d's with negative powers are zero and the 
only non trivial equation in the null space of M is that (i2,o = <^i,0) the c coefficients are, 

(2.44) 
(2.45) 
(2.46) 
(2.47) 
(2.48) 
(2.49) 
(2.50) 

Of course at this stage we in full agreement with the known results from Refs. [9, 19]. 

As a final remark we consider the possibility of fitting the integrand using the large 
momentum limit, t — )• cxd, matching to the method of Forde [10]. This derivation follows 
the argument presented in the recent review article of Ellis, Kunszt, Melnikov and Zan- 
derighi [6]. The subtraction terms, coming from box contributions previously evaluated 
with quadruple cuts, can be written schematically as, 

5(fe) = ^^^^^iif^, (2.51) 

where ^/^^i^i^i.^j^k') is given in eq.(2.10) using bj^{Pj) in eq.(2.3). The integrand is formed 
from the two scalar products. 



2fc« -(^(P,-) = [{K\\P,\K\\t-''-^{K\\P,\K\^ , (2.52) 
2fc(2) . ^(P,.) = [-^K\\P,\K\\t + ^(i^JlP,.]^-^]) . (2.53) 



-9- 



Explicitly taking the limit t — ?• oo yields, 

lim5(^)(t) = (-ir^ci,,,,,3, + Ori). (2.54) 

Since we know, 

3 

ri,i,i,(k) - s^'\t) = ds,xe (2.55) 

a;=— 3 

we can define a new set of coefficients d'^^^ 

3 

lim (7;,.,.3(^^'^)) =E<-*" + ^(*"')' (2.56) 

t— i-CXD \ / ^ ' ' 

X = 

such that d'g q = dsfi + ^-^^Ci;j^j2i3i ™d d'^ .j, = dg^x for all x > 0. Inverting gives. 

Coo = difi = ^2,0 = ^ (4,0 + 4,o) ■ (2-57) 

enabling us to extract the non-spurious coefficient from the triple-cut alone as derived in 
refs. [6, 10]. 

No further subtleties arise for double (or even single) cuts so we won't reproduce any 
further results at one- loop. 

3 Integrand Representations of Two Loop Amplitudes 

The extension of the generalised unitarity algorithm from one to two loops is complicated 
by the fact that a general loop integral basis is not known. However, it has be known for 
some time how to reduce a general two-loop Feynman diagram to a basis of master (no 
longer simply scalar) integrals by the use of integration by parts identities. For a cut based 
construction of the amplitude, such a basis is unfortunately not suitable since the doubled 
and crossed propagators that appear do not factorise onto simple poles, and hence neither 
onto the products of tree- level amplitudes. 

A unitarity compatible basis has been explored using a Grobner basis construction 
in a recent paper by Gluza, Kajda and Kosower [48]. Here will follow a slightly different 
approach in which instead of trying to fit the coefficients of a minimal basis of functions, we 
will fit an integrand level expression compatible with unitarity cuts, which can be further 
reduced to master integrals by any of the known techniques. This is similar to the integrand 
reduction technique recently presented by Mastrolia and Ossola [51]. 

For the purposes of this initial study we focus on the parts of the amplitude sensitive 
to seven propagator cuts in four dimensions. Though a small step towards a full integrand 
level reduction technique we will emphasise some of the features that we hope apply to a 
wider class of cuts. For the present paper we will be concerned with primitive amplitudes 
contributing to gluon-gluon scattering two loops. In particular this restricts us to the case 
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where no subtraction terms from octa-cuts is required though the procedure is expected to 
follow in a similar fashion [51]. The colour ordered partial amplitudes are defined by [33], 



^f) = (^6 tr(TMi)r'^.(2)T'^.(3)rM4))(^iV242),iC^^^^^^^^2),a(3),a(4)) 

cr€S4/Z4, 

+ (a(l),a(2),a(3),a(4))) 
+ 9s Yl ^ctr(rMi)rM2))tr(rM3)rM4))Ag^3((j(l),(j(2);a(3),(T(4)) (3.1) 



where Nc is the number of colours and T""^ are the fundamental generators of SU{Nc). 
These partial amplitudes are mapped to primitive amplitudes before unitarity cuts can be 
applied using colour ordered tree-level amplitudes. 



4}f{l, 2, 3, 4) = (1, 2; 3, 4; ) + Af""^^ (2, 3; 4, 1; ) + (1; 2, 3, 4; ; 

+ (2; 3, 4, 1; ) + 4^^°^^ (3; 4, 1, 2; ) + (4; 1,2,3;) 

4^J;f (1, 2, 3, 4) = 24^^^°^] (1, 2; 3, 4; ) + 24^^^°^! (1, 2; 4, 3; ) - 44<^^°^] (1, 3; 2, 4; ) 
+ 24'^"°^] (2, 3; 4, 1; ) - 44'^"°^! (2, 4; 3, 1; ) + 24^"°^^ (3, 2; 4, 1; ) 
+ 24"''"k3; 2, 1; 4) - 4Af ' (2; 3, 1; 4) + 24"''°"'(2; 4, 1; 3) 

+ 2Af "'(1;2,3;4) -4Af "^(1;2,4;3) +24"''°"'(1;3,4;2) (3.3) 
Ag,3(l, 2; 3, 4) = 6Af (1, 2; 3, 4; ) + G^"""^' (1, 2; 4, 3; ) + 24^'^°^! (1; 2, 3, 4; ) 
+ 24^"°"! (1; 3, 4, 2; ) + 2Af (1; 4, 2,3;)+ 2ylf '"^l (2; 1, 3, 4; ) 
+ 24^"°"! (2; 3, 4, 1; ) + 2Af (2; 4, 1, 3; ) + 24^^°^^! (3; 1,2,4;) 
+ 2ylf (3; 2, 4, 1; ) + 2Af (3; 4, 1, 2; ) + 24^"°^! (4; 1,2,3;) 
+ 2ylf (4; 2,3,1;)+ 2Af (4; 3,1,2;)+ 44^"°^1 (3; 2, 1; 4) 



. [dbox 



[dbox] 



1 [pbox] 



(3.2) 



24"''"'(2;3,1;4) -2Af'^°"](2;4,l;3) -24"''"'(1;2,3;4) 



2^[xbox] 2^ ^. 3^ ^ ^^ixooxj 3^ 4. 2) 



[xbox] / -, , 



(3.4) 



The decomposition involves three topologies: the double box, the crossed box and the 
penta-box which we reference explicitly in the superscript for clarity. We will examine 
the four dimensional hepta-cut part of these primitives in the following section. We label 
each integrand with a subscript for the number of cut propagators and a set of indices 
labelling the momenta leaving the diagram at each vertex. A '*' label indicates that no 
external momentum enters a vertex. Using this notation, which is described in more detail 
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in appendix D, the primitive amplitudes are written: 



A, iWA,) - j j (3.5) 

^ ^ ' ^"7 y (2vr)^(2vr)'^lv OLi ^1 J 

where ' . . . ' represents terms with < 6 propagators and terms only accessible via d-dimensional 
cuts. The above decomposition only applies to the pure gluonic loops but we may also use 
it for gluino and adjoint scalar loops. 

4 Hepta-cuts of Two-Loop Amplitudes 

We will proceed with the integrand reduction through a three step process which utilises: 

• Relations from Gram matrices to find a general form for the integrand. 

• Finding the total number of on-shell solutions, these will be families of solutions 
depending on a number of free parameters. 

• After fitting the full integrand, further reduction of non-spurious terms to Master 
Integrals (Mis) can be achieved using IBP relations 

The full integrand can then be constructed as the solution to a linear system of equations. 
In the following we go through the details of the three independent seven propagator 
topologies for four-point amplitudes with massless legs: the double box (shown in fig. 2), 
the crossed box (shown in fig. 3) and the penta-box (shown in fig. 4). 

4.1 Integrand Parameter isat ions from Gram Matrix Constraints 

Gram matrices are also important for two loop amplitude computation. Let /c, q be the 
loop momenta and {61,62,63,64} be the basis of 4-dimensional momenta. When k and q 
are 4-dimensional momenta, we obtain three quadratic relations of the Lorentz invariants 
using eq. (2.16), 

k'^ = {k ■ ei,k ■ e2,k ■ e3,k ■ e4)G^^{k • ei, /c • 62, A: • 63, A: • 64)"^ (4.1) 
q^ = (g • 61, g • 62, g • 63, g • 64)Gj^(g • 61, g • 62, g • 63, g • 64)'^ (4.2) 
p ■ q = {k ■ ei,k ■ e2,k ■ e3,k ■ e4)G'^^{q ■ ei,q ■ e2,q ■ e3,q ■ 64)^ (4.3) 
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Alternatively, when k and q are 4-dimensional-momenta, by the linear dependence property, 

A...detG^^^^^^^^^Uo (4.4) 

y ei 62 63 64 k J 

X,,^detG A = (4.5) 

\^ ei 62 63 64 g y 

A., = detG^^^^^^^^'Uo (4.6) 

y 61 62 63 64 q J 

It is easy to see that (4.4), (4.5) and (4.6) are equivalent to (4.1), (4.2) and (4.3), respec- 
tively. 

It seems that there are many more 5x5 Gram matrix relations. Because we can choose 
5 vectors from the set {61, 62, 63, 64, k, q} twice, there are 

l((^)+l)(^)=21 (4.7) 

5x5 Gram-matrix relations for 4-dimensional loop momenta. However, the additional 
relations are not independent, since they can be generated by (4.4), (4.5) and (4.6), by the 
linear and anti-symmetric properties of Gram matrices, (2.13) and (2.14). 
To see this, we consider k and q as general d-dimensional vectors, 

4 4 
k = "^ki6i + k'', q = ^qi6i + q'', (4.8) 
1=1 1=1 

where A:" and q"^ are the extra-dimensional components. It is clear that 

Xkk = detG ( "2 ! J = det(G4)(fc")^ (4.9) 
\^ 61 62 63 64 fc" J 

A,, = detG ( ^ J = det(G4)(g")2. (4.10) 

y 61 62 63 64 Q'" J 

Xkg = detG ( "2 ^ J = det(G4)(A;" • g"). (4.11) 

\^ 61 62 63 64 J 

For other 5x5 Gram matrices, comparing with the above expressions, 
idetci'' ^ A =k4detG('' f"V94detGf^^ 

y 61 62 63 64 g y y 61 62 63 64 A;" y y 61 62 63 64 fc" j 

+ det G ( ^2 '3 ? J = A:4Afc, - (?4Afcfc, (4.12) 
\ 61 62 63 64 A;" / 



detG ( ^1 ^2 63 A; \ ^ ^2_)^^^ _^ ^2^ _ ^i^^^^y^ _^ f ^^kk\q - AL)> (4-13) 

6\ 62 63 k q f det(G4'"^ ^ 



with other Gram matrices following a similar pattern. So as long as the relations (4.4), (4.5) 
and (4.6) hold, all the other 5x5 Gram matrices will vanish automatically. However, in 
practice, we will still use other 5x5 Gram matrix relations, since they usually provide very 
efficient ways of combining the three fundamental relations (4.4), (4.5) and (4.6) together. 
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4.2 The Double Box 




Figure 2. Conventions for ttie momentum flow in tlie planar double box. 



We define the double box contribution to A 



planar, (2) 



7dbox 
'1 

idbox 
'5 



-Q+P4 



;dbox 
'2 

;dbox 



k-pi 

-q 



dbox 



''3 

rdbox 



by the following propagators: 

k-pi,2 lf"''' = -q+P3A 
-q-k (4.14) 



Just as in the one- loop box topology, the three external momenta {^1,^25^4} are supple- 
mented by the spurious vector, u, given in eq. (2.3). The loop momentum is then contained 
in the space spanned hy v = {pi,P2,P'i,^}- Taking into account scalar products that can 
be trivially rewritten in terms of the propagators (^^''°^)^ using relations of the form of 
eq.(2.8), we can parameterise the integrand, Aj^°^^^^^^{k,q), with four irreducible scalar 
products (ISPs) combined into terms of the form. 



ik-p,r{q-p,nk.u;nq-u;f. 



(4.15) 



The first constraints on the indices m, n, a and /3 come from renormalization conditions 
implying m + n + a + /3 < 6. Since each of the integrals involves four propagators for this 
topology we may also deduce m + a < 4 and n + /? < 4. 

The Gram matrix relations (4.4), (4.5) and (4.6) put constraints on the ISP's. For the 
basis {61,62,63,64} = {pi,P2,P4,i^}, at the hepta-cut (4.4) reads. 



{k-ujf = {k-pi- su/2f 
This relation requires that a = 0, 1. Similar, (4.5) reads, 

{q-uf = {q-pi-si^/2f 
So /3 = 0, 1. (4.6) requires that a/3 = 0, since it reads, 



{k ■ (jj){q ■ (jj) 



§14 Sl4(,fe 

4 2 



Pi) ^ sii{q-pi) _^ A _^ 2si4 

S12 



{k-pi){q-pi] 



(4.16) 



(4.17) 



(4.18) 



The number of the ISP monomials is reduced to 56. 
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So far, we just used the three fundamental 5x5 Gram-matrix relation individually. 
It is possible to combine them together to get more constraints. The efficient way is to 
consider other 5x5 Gram-matrix relations directly. We have, 

detG =0, detG , = 0, detG ^ = 4.19 

\\2^kqj ' \\2^ujkj ' \\ 2 4 uj q j ^ ' 

For example, the first equation in (4.19) explicitly reads, 

= 4(/c-p4)^(g-pi)^ + 2si2(fc-p4)^(9 •pi) + 2si2(A;-p4)(g-Pi)^-si2Si4(A:-p4) (g-pi) (4.20) 

so the terms with both m > 2 and n > 2 are reduced. These relations further reduce the 
number of ISP monomials to 32. With all our constraints complete we arrive at a general 
parameterisation for the double box integrand, 

Afi?.34*(A;,'?)= (^mnia+2p){k-p,r{q-pink-unq-uf. (4.21) 

mnajS 

There are 16 non-spurious terms, i.e. those not proportional to {k ■ to) or (q ■ uj), 

(COOO) CoiO, ClOO) C020) Clio, C20O) C030, C120, C210, C300 , Co40 5 C130 , C310 , C400 , C140, C410) (4.22) 

and 16 spurious terms, 

(cooi) Coil, cioi. Cm, C201, C211, C301, C311, C002, C012, C102, C022, C112, co32, C122, C132) (4.23) 
The terms can be represented in form of a table: 

a = 0, /3 = 





n=0 


71=1 


n=2 


n=3 


n=4 


m=0 


/ 


/ 


/ 


/ 


/ 


m=l 


/ 


/ 


/ 


/ 


/ 


m=2 


/ 


/ 








m=3 


/ 


/ 








m=4: 


/ 


/ 









and for the spurious terms, 

a = l, /3 = a = o, /3 = 1 





n=0 


n=l n=2 n=3 




n=0 


77 = 1 


n=2 


71=3 


m=0 


/ 


/ 


771=0 


/ 


/ 


/ 


/ 


m=l 


/ 


/ 


m=l 


/ 


/ 


/ 


/ 


m=2 


/ 


/ 


m=2 










m=3 


/ 


/ 


771=3 











Our next task is to use the full set of on-shell solutions to find a map to these coefficients 
from the products of tree-level amplitudes 
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4.2.1 Solutions to the on-shell constraints 

The solutions to the on-shell constraints = have been considered in Refs.[50, 51]. 

The six solutions can be parameterised using the same two component Weyl spinor basis 
as used at one- loop: 

^2 = ^iPl + ^•2P2 + XZ ^ h Xi 

I5 = ym + y^Pi + ys — ^ — + y4 — - — (4.24) 

With eight unknowns and seven equations, each of the solutions depends on a free param- 
eter which we will denote as r. The choice of this parameter has been made to ensure that 
the integrand takes a simple polynomial form. 



Solution 


Xi 


X2 


X3 


X,i 


yi 


y2 


ys 


yi 


1 
2 
3 
4 
5 
6 


















(23> 
(13) 



[24] '^^ ^) 




(13>(24>r+(12>(34> 
(14)(13)r 




[23] 
[13] 



(24) V 
[13][24]r+[12][34] 
[14][13]r 




















(13) ^) 



[23] 
[24] 






(23) 
(24) 





It is straightforward to feed these solutions in the general integrand expression A7;i2*34*(A;, q) 
and define a set of coefficients that can be extracted from the product of six tree level am- 
plitudes, 

A flis) (s)\ a(«) t \ )Y!,t=0'^f!,xT^ 1,2,3,4, rAOK\ 

A7;12*34*(A;^ ',q'' ')= A^.{2*34*r) = ) , ^ . „ (^-^^j 

where, 

A7;i2*34*(g,A:) = A^''\-l^^\p,,l^')A^^\-q^\p2,l^nA^'\-q^\p3,l^') 

Afc=± 

X AW(-/5-^^P4,/6'^)^(°^(-/6^^^/^,/7')^^°^(-^3"^^^4^-^7"^')• (4-26) 

We now follow the procedure used in section 2.2 by constructing a 38 x 32 matrix such 
that, 

d = M -c (4.27) 

It is easy to check that this matrix has rank 32 and therefore a unique solution. We are able 
to invert the system using standard linear algebra packages available for symbolic compu- 
tations. The final list of equations mapping ds^x to c^.n{a+2(5) is available in a computer 
readable format from http://www.nbia.dk/badger.htnil. They have relatively simple 
forms for example: 

cooo = 2 (^1,0 + d2,o) ■ (4.28) 
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4.2.2 Integration by parts identities 

Having obtained a method to fix the 32 coefficients of Aj-i2*34*{k,q), it is now in a form 
that can be further reduced to master integrals using integration by parts identities. There 
are by now a number of packages available to perform the task of reducing the master 
integrand, ■ 12*34,* {k, q), onto a basis of two master integrals. For this purpose we made 
use of the FIRE Mathematica package [57]. In the case of the planar double box this enables 
to compare our results directly with those of Kosower and Larsen [50] . 



4'^^°^1(1,2;3,4;)= J J 



d'^fc d'^g Ci + C2{k-P4) 



+ 



(4.29) 



where we suppress all further master integrals. In terms of our non spurious CnmO coeffi- 
cients they are, 



Ci 



, S12S14 

COOO H n Clio 



2 3 

(Cl20 + C210) H 5^ (Cl30 + C310) 



16 



32 



64 



C2 — ClOO + ColO 



(Ci40 + C410) , 
3S12 



3^2 

Clio + -ir (C020 + C200) + — ^ (Cl20 + C210) 



£14 

2 



+ 



'14 



(C030 + C300) 



— TT (Cl30 + C310) + 

16 



3 
14 



,C040 + C400J 



+ (Cl40 + C410J 



(4.30) 



(4.31) 



Closed form expressions for the master integrals can be found in refs. [61, 62]. 
4.3 The Crossed Box 




Figure 3. Conventions for the momentum flow in the non-planar crossed box. 



We represent crossed box topology shown in fig. 3 using the propagators as follows, 

k 



/xbox 
'1 

/xbox 
'5 



k+pi 

q-P4 



/xbox 



/xbox „ 

'4 —I 



xbox 



k + P2,3 



1^ 
«7 



xbox 



q-k+p2, 



(4.32) 
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As for the double box there are four ISPs parameterising the general integrand, A7;i*34^,2(A:, q), 
the Gram matrix constraints however lead us to a slightly different form of the final in- 
tegrand. Again, (4.4), (4.5) and (4.6) confine (a,/?) = (0, 0), (1, 0), (0, 1), respectively, 
giving, 

A^,Z4Ak,q) = Yl ^mnia+2(S)ik-P3riQ-P2nk-u;nq-ojf. (4.33) 

mnajS 

We can combine the three fundamental relations together or consider other 5x5 Gram 
matrix constraints, to reduce the dependent terms. These relations have different form 
comparing with the double-box case, since there are fewer symmetries in the crossed-box 
diagram. This leads us to a representation with 19 non-spurious terms, 

(COOO ) COlO ) Cioo , C020 ) Clio, C200) CqSOj C1205 C210 , C300 , Co40 j C22O j C310, C400, C050, C320, C410, CoeOi C420) 

(4.34) 

and 19 spurious terms, 

(C001> Coil, Cioi, Cm, C20I, C211, C301, C311, C002, C012, C102, C022, C112, C032, C122, C042, Ci32, C052, C142) 

(4.35) 

The terms can be represented in form of a table: 



a = 0, 13 = 





n=0 


n=l 


n=2 n=3 n=4: n=5 n=6 


m=0 


/ 


/ 


/ / / / / 


m=l 


/ 


/ 


/ 


m=2 


/ 


/ 


/ 


m=3 


/ 


/ 


/ 


m=4 


/ 


/ 


/ 



and for the spurious terms. 







a = 


1, /3 


= 






n=0 


n=l 


n=2 


n=3 n=4 


n=5 


771=0 


/ 


/ 








m=l 


/ 


/ 








m=2 


/ 


/ 








m=3 


/ 


/ 












a = 


0, (3 


= 1 






n=0 


n=l 


n=2 


n=3 n=4 


n=5 


m=0 


/ 


/ 


/ 


/ / 


/ 


m=l 


/ 


/ 


/ 


/ / 




m=2 












m=3 
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4.3.1 Solutions to the on-shell constraints 

The crossed box topology uses the same basis as the planar box, 

q = xip'l + X2P'^ + X3 h X4 

^4 = yiPs + y2P4 + 2/3 ^ + 2/4 ^ 



(4.36) 



which we then use to for solve {/^} = 0. 

The result is 8 families of solutions, again parameterised by r. These can be sum- 
marised by: 



Solution 



1 

2 
3 
4 
5 
6 
7 



Xi 


3^2 


S14+T 





Sl2 


Sl4+r 





Sl2 














S14+T 





Sl2 


S14+T 





Sl2 


-1 





-1 






(23)(si3+r) 
(13)si2 






[13](£i4+t) 



[23]si2 



Jjl4+r 
(13) [32] 



X4 





[23](si3+t) 

1131^12 





[13] (32) 
(13)(si4+t) 



(23)si2 



SH+T 

[13] (32) 





yi 



2/2 



2/3 



P2yi24i 



_[23] 
[24] 





[24] (32) 



[32] 
[14] 





2/4 









_M 

(24) 





(24) [32] 



M 

(14) 



From which, thanks to the choice of r in eq. (??), we define the polynomial form of the 
cut integrand. 



A7;i*34*2(A:(''),g('')) = A 



7:1*34*2 



{r) 



1,2,5,6, 



E'=o^^,--^" 5 = 3,4,7,8. 



(4.37) 



where, 

A7;l*34*2(g,fc) = A^'H-ll'\Pl,l2nA^'H-le'^P2,l7nA^''\-l3'\P3,lA') 
Afc=± 

xA(o)(-/4-^^P4,^5'')^(°H-/5"'^^l^^6'')^^°^-^2"'^^3^-^7"'o■ (4-38) 

As before this leads to an invertible matrix this time 48 x 38, 

d = M c (4.39) 

The final equations for Cnm{a+2p) of similar complexity to those in the double box 
topology for example: 



cooo — 2 (^1.0 + "^2,0) • 



(4.40) 



The complete set of relations can be obtained from http://www.nbia.dk/badger.html. 
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4.3.2 Integration by parts identities 

The integration by parts identities generated using FIRE reduce A7;i^,34*2 onto two seven 
propagator master integrals, 

where Ci and C2 are given by: 

Cl = COOO + T7:'S14S13(C200 — Clio + 2co2o) 
lb 

+ — S14S13(S14 - S13)(C300 " C210 + C12O " 2co3o) 

+ Y^(3(S14 - Svif + Si2)si4Si3(c400 - C310 + C22O + 2co4o) 

+ Y^(('^14 - Sl3)^ + Sl2)'5l4Sl3(si4 " Sl3)(c320 " C410 " 2co5o) 

+ Y^(5(S14 - Svi)'^ + 10Si2(si4 - Si3)^ + s\2)SiiSvi{Ci20 + 2co6o) (4.42) 
3 

C2 = ClOO - 2C010 + -(Sl4 - S13)(C200 " cno + 2co2o) 
o 

+ :i77(2(si4 - Si3)^ + Si2)(c300 " C210 + C12O " 2co3o) 

lb 
2 

+ Y^(5(S14 - Sl3)^ + 7^12) (Sl4 - S13)(C400 " C310 + C220 + 2co4o) 
+ Y^(3(S14 - 513)"^ + 8Si2(si4 - Sl3)^ + 42)(c320 " C410 " 2co5o) 

+ (7(S14 - S13)^ + 30S?2(S14 - S13)^ + '^^si^) (su - Si3)(c420 + 2co6o) (4.43) 

The integrals themselves have been computed using Mellin-Barnes techniques in refs. [63, 
64]. 

4.4 The Penta-Box 

Our conventions for the penta-box topology follow those outlined in fig. 4. The seven 
propagators are, 

li = q - k-pi l2 = q-k+ p2,3 ^3 = 1? + P2,3 h = Q + P3 

k = q k = q-P4 h = k (4.44) 

This topology has a rather different integrand structure than our previous cases. The Gram 
matrix relation, 

(fc.u;) = 2(g.a;)f^^-^^') (4.45) 

V ■512 Sl4 / 

implies that [k ■ uj) is not independent at the hepta-cut. So we only have three ISPs, 
(^ ■ P2)) {k ■ P4) and (q ■ oj), with a final form parameterisation, 

A?S234* = E ^--(^ • • P^riQ ■ (4.46) 

m,n,a 
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4 1 




Figure 4. Conventions for tlie momentum flow in the penta-box. 
The sums are restricted such that there are 20 coefficients in all. Ten are non-spurious, 

{COOO) ClOO, CoiO, C020> Clio, C200; C030, Cuo, C210> C300}, (4.47) 

and ten are spurious, 

{cooi) cioi, Coil, C021, Cm, C201, cosi, C121, C211, C301}. (4.48) 
In the tabular format this looks like, 

a = a = 1 





n=0 


n=l 


n=2 n=3 




n=0 


n=l 


n=2 n=3 


m=0 


/ 


/ 


/ / 


m=0 


/ 


/ 


/ / 


m=l 


/ 


/ 


/ 


m=l 


/ 


/ 


/ 


m=2 


/ 


/ 




m=2 


/ 


/ 




m=3 


/ 






m=3 


/ 







4.4.1 Solutions to the on-shell constraints 

We parameterise the loop momenta according to, 

q = Xip'l + X2P'2 + X3 h X4 

I5 = ym + y-^Pi + 2/3 — ^ — + y4 — ^ — • (4.49) 

Interestingly in this case we find that the set of cut constraints is degenerate and we have 
two independent solutions parameterised by ri and T2, 



Solution 


Xi 


X2 




X4 


yi 


y2 


ys 


yi 


1 
2 


Tl 






(23)(S12S14+T1+T2) 




[23](si2Sl4+Tl+T2) 










[23\ 
[24] 






(23) 
724Y 


S12S14 
Tl 


(13)si2S14 




S12SH 


|13|si2Sl4 
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The choice of ri and T2 has been made such that the integrand has a symmetric form for 
the ISP's and we are able to write down a simple form for the inverted 20 x 20 system: 

Cmna = (_ 1)™ 2"+-+"-l ^MIk{di,„,,n + (-l)"ci2,m,n) (4.50) 

where ds^m,n is the coefficient of T{"r2 for solution s. 
4.4.2 Integration by parts identities 

After the application of further reduction via IBP relations it turns out that all penta-box 
integrands are reducible to six propagator master integrals or simpler topologies. Since 
those master integrals will also have contributions from hexa-cut configurations a com- 
plete study of the form of these reduction identities will be postponed to future studies. 
Nevertheless, for the purposes of a complete integrand level reduction these terms play an 
essential role. 



5 Applications to Gluon-Gluon Scattering 

In this section we apply our technique above to gg — )• gg scattering amplitudes. These am- 
plitudes have been known from some time [40, 41] and present something of a benchmark 
test of our method. We compute the amplitudes in Yang-Mills theory with an arbitrary 
number of massless gluinos (nj) and scalars (n^) in the adjoint representation. By consid- 
ering specific configurations of the number of fermion and scalar flavours we are able to 
cross check our results against the simpler ones obtained in super-symmetric Yang-Mills 
theories, using: 

(nj = 4,n, = 3)^AA = 4SYM, (5.1) 
(?!/ = 2,n, = 1) = 2 SYM, (5.2) 

(nj = l,n, = 0) ^AA= 1 SYM. (5.3) 

Note that we consider scalars to be complex so there are two degrees of freedom for each 
scalar flavour. 

Throughout this section we will use 

A?:,....,(A:,g) = (Ag:.,^+A?:t..,JAW, (5.4) 

where Ap^^,^ contain all spurious dot products, which vanish after integration, and 
■^T-iT 16 contains the remaining non-spurious dot products. The four-point tree ampli- 
tudes shorthand for the functions A^^\l^^ , 2^^ , 3^'^ , A^*) 

given in Appendix B. 

There are ten configurations of particles flowing in the loops which contribute to the various 
topologies, the explicit cases of the double box are shown in fig. 5. 

The final results are factorised into terms that vanish depending on the amount of 
super-symmetry. We note that in = 4 the factors (4 — nj) , (3 — n^) , (1 — -|- n^) 
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4 14 14 1 




Figure 5. Contributions to the double box in Yang- Mills theories of gluons, Uf adjoint fermions 
and ns scalars. Gluons are represented as solid lines, fermions as solid lines with arrows and scalars 
as dashed lines. 

all vanish whereas (l — nf + Us) will only appear in theories with no super-symmetry. 
Formulae as a function of the number of super-symmetries, M, can be obtained by: 

nf=J\f, ns=M-l. (5.5) 

For the planar and non-planar double boxes we have cross checked the final coefficients 
of the master integrals are in full agreement with the results by Bern, De Freitas and Dixon 
[41]. 

5.1 Planar Double Boxes 

The following section contain the results for the planar double box as defined in eq. (4.21). 
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5.1.1 The h+ Helicity Amplitude 

We find the fohowing analytic forms for the integrands: 

Adbox,ns, h+ 2 /r r\ 

7;12*34; = -S12SU (5.6) 

A^S;-"-^ = (5.7) 

After applying further IBP relations as given in eq. (4.31), the coefficients of the two basis 
integrals [48, 50] become: 

Cr-++ = -4^14^^°^ (5.8) 
C,-++ = (5.9) 

5.1.2 The - H h Helicity Amplitude 

The integrands for this configuration are: 

Adbox,ns, — I h 2 
7;12*34* - -Sl4'Sl2 

(4 - Uf) (3 - Hs) Si4Si2^ 



{{k ■4) + {q- 1)) {2{k ■ 4) {2{q ■ 1) + su) + su {2{q ■ 1) - sm)) 



Sl3 ^ 

+ 4{k ■ 4) (6si2(g • if + S12 (2si2 - S14) (<?•!)+ 8((? • 1)^ - si2Si4^) 
- 16{k ■ 4)4 - (si4 (2((/ • 1) + S12) - 4(g • l)^) 2' 

(4 - Uf) gl 4Sl2^ f^^,^ ^^2 

+ 4 (3si2 + S14) [q ■ if - 2si4 (3si2 + S14) (? • 1) + S12S13S14) (5.10) 

A dboXjS, — \ h 

^7;12*34* — 

^^^"""^^^^ "3"^^ ''''''' • 1) • 4) + ^12) (g • -) - (A. • 4)(fc • c) (2(, . 1) + .12)) 



+ ^ . 3 4(fc • 4)2 (4(g • 1) + 3si2 + 514) + 2{k ■ 4) {2{q • 1) - su) (4(g • 1) + 3^12 + su) 



+ 



Sis" 

514512^ jl-nf + Us 



(A; • io){^{k ■ 4)2 (4(g • 1) + S14) + 8(fc • 4) (3si2(g • 1) + su (2(g • 1) + S12)) 

- 16(A; -4)3 + 512' (2(^7 •l) + si4)) 

+ {q-uj){-2{k- 4) (4 (3si2 + 2si4) {q-l) + 16(g • 1)^ + S12') 

- S14 (8si2((/ • 1) + S{q ■ if + S12') + 16(g • 1)3)) 

+ ^^^7^^(«12513(-(A^ • CO)) {2{q ■ 1) + su) 

+ 2{k ■ 4) (2si4(A: • to) (4(g • 1) + 3si2 + su) + (s^sis - 8su{q ■ 1)) {q ■ w)) 



+ si4 (si2Si3 - 4 (3si2 + S14) ((? • 1)) {q-u})) (5.11) 
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which lead to coefficients of the master integrals of, 



^ W^^(o) (2,^^,^^ il-nf + Us) - (4 - n;) .13') (5.13) 



5.1.3 The - + H — Helicity Amplitude 

The integrands for this configuration are: 



Adbox,ns, — \ — I — 2 
7;12*34* - -•S14S12 

- - + ((k . 4)3 (8(<7 • 1) - 4si4) + (A; • 4)2 (.14^ - 2 (3.12 + 4.14) {q ■ 1)) 

+ (A; • 4)(g • 1) (4((7 • 1) - 3si2 - 2^14) {2{q ■ 1) - su) + 4(/c • 4)^ + (g • 1)^ (5^4 - 2{q ■ 1)) ^ 
+ ^(i^M(^p£)fiH!(, . 4)(, . 1) (2((. . 4) + (, . D) - .14) 

+ ^^^^4^(si4'((A: • 4) + (<? • 1)) - 2.14 {-{k ■ 4)(g • 1) + (A: • 4f + {q • l)^) 

Sl4^ V 

- 8{k ■ A){q ■ l){{k ■ 4) + {q ■ 1))) (5.14) 



A dbox,s, — \—\ — _ 
^7;12*34* — 

(4_^ )(3_^^)g^,2 / . ^) + . a;)) - 2si2.i4((A: ■ 4){q ■ u) + {q ■ l){k ■ u)) 

2si3Si4^ V 

+ 8si3{k ■ 4){q ■ l){{k ■ oj) + {q ■ u))^ 

+ ^"'^^"V"^^ (^W((A:-u;) + (,..)) 

+ si2Si4^ (3si2((A; ■Lo) + {q- u)) - 4{k • 4)(g . u) - 4(g . l){k ■ uj)) 

+ 2si4(4si3 {2{k . 4)(g . l){{k ■ lo) + {q ■ oj)) + {k ■ i)\k ■ lo) + {q ■ l)\q ■ oj)) 

- 3si2\{k ■ 4)(g ■uj) + {q- l){k ■ to))) - Ssrsi " ^^i^ik ■ 4)(g . l)((fc ■ u) + (q ■ co)) 

+ 2{ik. ifiik ■ 4) + 2{q . l)){k .u;) + iq. l)\2{k ■ 4) + (q ■ l)){q -Oj)))) 

+ ( _ 3si2Si4'((A: ■io) + {q- u)) - Ws^ik ■ A){q ■ l){{k ■ u) + {q ■ co)) 

2si3Si4^ V 

+ 2si4 (3si2((fc ■ 4){q ■ u) + (q ■ l){k ■ u)) - 2si3((A; ■ 4)ik ■ co) + (q ■ !)((/ • co))) ) (5.15) 
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which lead to coefficients of the master integrals of, 



2 
2_ 

2\ 



Ci-++- = (2si2 (I0si2^ + 11S14S12 + 2si4^) (1 - n/ + 



+ si4 ((4 - n/) (3 - n^) S12 (2si2 + su) - (4 - Ji/) S12 (4si2 + su) + 4si4 ) ) 

(5.16) 

C2-++- = ( (2OS12' + 22S14S12 + 4514^) (l-Uf + ns) 

ZS14 

+ si4 ((4 - nj) (3 - n,) (2si2 + sm) - (4 - n/) (4si2 + su)) ) (5.17) 

We notice that the explicit expressions for all helicity amplitudes never contain tensor 
coefficients higher than rank four. In other words the coefficients C410, C140, C311, and C132 
are zero even in pure Yang-Mills though we were not able to exclude them a priori from 
the renormalization constraints. 

5.2 Non-Planar Crossed Box 

5.2.1 The h+ Helicity Amplitude 

The integrands for this configuration are: 

Axbox,ns, 1 — h 2 
7;1*34*2 - -•S14S12 

+ i (4 - Uf) si4 (-2 (si2 + 2si4) (g • 2) - A{q ■ 2f + S13S14) 

+ '^'^^-''1 + ''^^ ( _ 4(5 . 2f {2{q . 2) + su) ' - ^14^ {-8s,s{q ■ 2) + 24(g • 2)^ + s,^) 
+ 4sii{q ■ 2) (-6si2((? • 2) - 8{q ■ 2f + S12S13) ) (5.18) 



A xbox,s, 1 — h 

^7;1*34*2 — 

- ~ ((g . u) {si3 {2{q ■ 2) + S14) - 2si2(fc • 3)) + 2{k ■ w) ((si2 + 2si4) (<? • 2) - S13S14)) 

+ (1 - rif + ?!,) / _ (5^2'(3(A; • w) - (<? • w)) + 3si4Si2(2(A: ■ u) + {q ■ w)) + 6si4'((7 • oj)) 
+ 16 (si2 + 2si4) (g • 2f{q ■ uj) + lQ{q ■ 2f{q ■ w) + si2Si3Si4((g • " 2(A: • w)) 
+ 2(A; • 3) (16 (si2 + 2si4) {q-2) + 16(g • 2)^ + 5^2^ + 12si4^ + I2S12S14) {q ■ c^)) (5.19) 

After applying the IBP relations we obtain the following coefficients of the master integrals 
via eq. (4.43), 

Cr^- = ^.14^(0) ^_ W£M^il_n^ ^ _ ^^^^^^ _ 4^^^2^ (5.20) 
^2""''"' = -ip^A^^^ (si3 - 514) ((4 - nf)si2^ - 2si3Si4(l -nf + n,)) (5.21) 
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5.2.2 The - H h Helicity Amplitude 

The integrands for this configuration are: 



Axbox,ns, — I h 2 
7;1*34*2 - -•S14S12 

^ l-i — ibf } CJT/IOT-^ 



(±J^f}^^f}lf _ 2.,3 (3(fc . 3){q • 2) + 3(^^ • 3)^ + {q • 2)^) 



+ 8(A; • 3){q • 2){{k ■ 3) + (q • 2)) + sis^q • 2) 
2(4-n/) (3 - n,) si4Si2^ 



(A;-3)((A;-3) + ((7-2))(2(g-2)-si3) 



4si4Si2^(l-n/ + 



n. 



+ ^ 4 ' 2su{k ■ 3){{k ■3) + {q- 2)) {2{q ■ 2) - S13) - 8(A: • 3)^(g • 2) 

Sl3 ^ 

- 4ik ■ 3f{q ■ 2)2 - 4(fe • 3)4 - si^^q ■ 2f + 4^13(9 • 2)^ - 4(g • 2)^) (5.22) 



A xbox,s, — I h 

'^7;1*34*2 — 

(4-n;)gi2gi4 / ^ _ (3si3(2(A: • 3) + • 2)) - 8(A: • 3)((/ • 2)) 
2si3-^ V 

+ ((? • w) (2S12S13 (-3(A; • 3) - 5(g • 2) + S13) + S14 (8(A; • 3) + 3^13) (si3 - 2(g • 2))) 
(4 - Uf) (3 - Us) S12S14 



^ 3 ,2si2(A;-w)(si3(2(A:-3) + (g-2))-4(/c-3)(<7-2)) 

+ {q-uj) {2{k ■ 3) (2si4 (si3 - 2{q ■ 2)) - suSis) + 513^ (2(g • 2) + su 



si4 (1 - n/ + 



n. 



+ (25i2(fc • ^)(85i2(fe • 3)2((A; • 3) + ((/ • 2)) 

+ S14' (4(A: • 3)((7 • 2) - si3(2(A; • 3) + • 2))) 

+ 2si4 (si3 (si2(2(fc • 3) + (<7 • 2)) + 2(A: • 3)^) - 4{k ■ 3)(g • 2) ((fc • 3) + s^)) ) 

+ (g • oj){su^ (- (8(A: • 3) (si3 - 3{q ■ 2)) + 513^)) 

+ 2514^ {{k ■ 3) {8{q ■ 2)2 + 3S12S13) + Si3^ {si2 - {q ■ 2))) 

+ 4S12S13514 {sisiq ■ 2) - si2{k ■ 3)) - isu\q ■ 2) (si3 - 2{q ■ 2)) 2)) (5.23) 
which lead to coefficients of the master integrals of, 



Cr+-+ = (2512514^ (1 - + n,) - (4 - n^) S14S13' - 4^13^) (5.24) 

C2^~^ = fi^^W (^13 + 3si4) (2S12514 (1 - n; + n,) - (4 - n;) ^13') (5.25) 
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5.2.3 The - + H — Helicity Amplitude 

The integrands for this configuration are: 




a: 



xbox,s, — I — I — 
■7;1*34*2 



^- ^^l±^(2si2{k ■ ^)(5i4' m ■ 3){q ■ 2) - si3(2(A: ■ 3) + (q ■ 2))) 

+ 2si4 {si3 {si2{2{k ■3) + {q- 2)) + 2{k ■ 3f) - 4(A: • 3)(g • 2) ((A; • 3) + S12)) 

+ 8si2(/c • 3) {{k ■ 3){{k ■3) + {q- 2)) - s^iq ■ 2)) ) 

+ (q ■ oj){su^ {8{k ■ 3) {3{q ■ 2) - S13) - sig") 

+ 2514^ {{k ■ 3) {8{q ■ 2f + 3S12S13) + si2 (5si2(g • 2) + sis^)) 

+ 4si2'si4 {{q ■ 2) {-A{k ■ 3) + i{q ■ 2) + S12) + Si^ik ■ 3)) - 2si4^g • 2) + IGsu^q ■ 2)3)) 



= 7^^^°^ (2si2Si3'(l -uf + ris)- ((4 - 7i/)si3 + 4si4)) (5.28) 
^2^^^^ = (3^13 + S14) ((4 - n/)si42 - 2si2Si3(l -nf + n,)) (5.29) 



As in the planar double box, non-zero tensor coefficients never appear higher than rank 
four since the values of the coefficient C050, C320, C410, coeo; C420, C052, or C142 turn out be 
zero independent of the helicity configuration. 

5.3 Penta-Box 

This topology only appears in helicity configurations which are zero at tree-level. There 
are two independent contributions, both of which vanish in super-symmetric theories. 




(2si2(fc • u) (si4(2(A; ■ 3) + {q ■ 2)) + 4{k ■ 3){q ■ 2)) 




(5.27) 



which lead to coefficients of the master integrals of. 
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5.3.1 The- + - 



— Helicity Amplitude 



^tS'sS ^ =4f(l-nj + n,) 
(13)2(14)2s3 
712; 



X -^'^ (fc • 2)(A: • 4) - —(A: • 2)2(A: . A) + —{k ■ 2)ik ■ 4)^ ) (5.30) 

(12)^sf3Si4 V ■512 Su 



(13)'(14)2sf2 ^ . 2)(A; • 4) - —{k ■ 2f{k ■ 4) + —{k ■ 2){k ■ Af ] {q ■ u) (5.31) 



5.3.2 The \ — Helicity Amplitude 



AS:234r-"- = 4.(i-n,+ 

X ^T^S^^''"' {^^ • 2)(fc • 4) - ^{k . 2f{k •A) + ^{k-2){k- Af) (5.32) 

(13)^Si2Sl4 V ^12 Si4 J 

A7S234* = (1 - + ^s) 

/12\2/14\2„2 / 2 2 \ 

^ /iU (^^ • 2)(fc • 4) - —(A. • 2f{k . 4) + — (fc • 2)(fc • 4)2 (g • ^) (5.33) 

(13)2si2Sf4 V S12 Sl4 J 

6 Conclusions 

In recent years, unitarity methods have been particularly useful in the computation of 
multi-loop scattering amplitudes in super-symmetric gauge theories and gravity. In this 
paper we considered the possibility of computing two-loop scattering amplitudes in a gen- 
eral renormalizable gauge theory with no super-symmetry via generalised unitarity cuts. 

The traditional unitarity approach to one-loop amplitudes relies on knowing a basis of 
scalar integrals in advance of the computation. Since such a basis is not known at two-loops, 
we looked to Gram matrix identities to constrain the general form of the integrand. This 
polynomial form can then be efficiently fitted by systematically evaluating products of tree- 
level amplitudes over a complete set of complex on-shell solutions to the loop momentum 
cut constraints. We derived a general map between the expansion of the tree level input 
and the coefficients of the integrand using only elementary linear algebra. 

The general integrand can be reduced to a set of master integrals using well known 
integration by parts identities. Using such identities we have derived master formulae for 
the three independent seven propagator topologies for 2—7-2 scattering. The method 
applies equally well to planar and non-planar topologies. 

As a test of our approach we computed the hepta-cut part of two-loop helicity ampli- 
tudes in Yang- Mills theory with adjoint fermion and scalars. This allowed us to check our 
results against the known results in super-symmetric Yang-Mills theories. 

Though a small step towards the complete reduction of an arbitrary two-loop ampli- 
tude, we hope the Gram matrix method introduced here will be of use in studying both 
D-dimensional cuts and cuts with fewer propagators. The extension to treat amplitudes 
with a higher number of external legs should also be possible following the basic steps 
described here, nevertheless a large number of basic topologies would be required. 
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Another interesting direction would be the appHcation of the technique to higher loop 
amplitudes. Though the solution to the Gram constraints will certainly be much more 
involved, the basic procedure for parameterising the integrand would be expected to apply. 
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A Conventions for Spinors and Spinor Products 

Throughout this paper we will make use of the well known spinor-helicity formalism in 
four-dimensions. A massless vector, p, can be written as, 

p'' = \{p\r\p]- (A.i) 

where are {p\ and \p] two component Weyl spinors of negative and positive helicity respec- 
tively. Helicity amplitudes can then be written in terms of spinor products {pq) and [pq] 
where Spq = {p + q)"^ = {pq) [qp\- We have also made use of the massless decomposition of 
two massive vectors. Pi and P2, in a basis of two massless vectors, K\,K2, 

- 7?2 - P?Pi ^ ^ ^ 

j^,,, _ 712 (712P,^ - PlPf) 

-'^2 - ^,2 p2p2 y^-'^) 

7l2 ^2 



712 = Pi-P2 + sign(Pi • P2) V (Pi • ^2)2 - PfPl (A.4) 

where we choose the sign in front of the square root to ensure that 712 > 0. Using these 
definition the size of the spurious vector in eq. (2.3) is, 

2_ detG(Pi,P2,P3) 
-~ detG(Pi,P2) ■ ^^-^^ 
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B Tree Level Amplitudes 



For completeness we list the well known formula for the independent tree level helicity 
amplitudes used in this paper: 



^(°)(l-,2-,3+) 


_^ (12)3 
(23)(31) 


(B.l) 


^(°)(l-,2-,3±) 


_ (12)^ 
- (31) 


(B.2) 


2",3r) 


_^(23)^ 
(31) 


(B.3) 


^ (Is ) 2 , 3s) 


.(12)(23) 
" (31) 


(B.4) 


^ (Is ) 2 , 3s) 


.(12)(23) 

" (31) 


(B.5) 


^^*^^(lg > 2g , 3s) 


= i(12) 


(B.6) 




= i(12) 


(B.7) 



the MHVamplitudes are obtained by complex conjugation (() -H- — []). All other ampli- 
tudes not related by parity or cyclic symmetries are zero. 

The non-zero four-gluon amplitudes used in section 5 are: 



^(0)(l-,2- 


,3+, 4+) 


AW(1-,2+ 


,3+, 4-) 


A(0)(l-,2+ 


,3-, 4+) 



(12)^ 



(12)(23)(34)(41) 

(14)^ 
^(12)(23)(34)(41) 
(13)4 



(12)(23)(34)(41) 
C General Solution for the One-Loop Box 



(B.8) 
(B.9) 
(B.IO) 



Just the for the purposes of completeness we give the explicit expressions for the terms 
used in section 2.1. We take the basis {Pi, P2, P3,u)} to span our space as before. For 
arbitrary kinematics with massless propagators the solutions are: 

k^'^'' = o '202 (^1 (712 + Pi) - Pi (712 + Pi) K'f 

_ ^ + ^^ kAmK\\ - kK\mK\^ (C.l) 

k^'^'" = o 'pJ Pl (712 + Pi) - Pi (712 + Pi) kY 

" ^ ;{K\h^\K'2] - .l!:?.. (i^^l7li^t]) (C.2) 
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where the massless vectors K\^, K\ are defined in Appendix A. The constant C is 

C =2P| (712 + Pf) {K\ ■ P3) - 2Pl (712 + Pi) {Kl ■ P3) 

- [1I2 - PM) {Pi + 2(^'2 • P3)) (C.3) 

while D can be expressed as: 

D = C^- iPiPi{Kl\P3\K^^]{K^^\P3\K{] (712 + Pi) (712 + Pi) (C.4) 
which is related to the normalisation of the spurious coefficient C2 by, 

VV4 = - (7?2 - PiPi) (C.5) 
D Notation for the Two-Loop Integrands 




lr-1 




lr-1 



(a) 



(b) 



Figure 6. Notation for a general ordered two loop amplitude, (a) A topology with two intersec- 
tions and (b) A topology with a single intersection. 

An n-point two-loop primitive amplitude can be described by an ordered set of mo- 
menta on external lines, say S = {pi,-" ,Pm} and another set on internal lines, T = 

{Pm+l,--- ,Pn}- 

We label the integrands by a set of indices, {ix}, corresponding to the position in the 
set S and a set of indices, {jx}, for the position the set T. The momentum leaving each 
vertex x = 1, . . . ,t is given by: 



ix+i — i 



(D.l) 



where the sum is considered to be cyclic modulo m. The momentum leaving each vertex 
y = t + 1, . . . , s are given by: 



Py = E Ta, 

a.=jk 



(D.2) 
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where jt+i — 1 = n. 

The vertices with intersections between the loops require a shghtly more elaborate 
notation. Each intersection is labeled by a set of indices give the ranges of momenta 
entering at each section, [ix,jy]- 

The integrand functions which we call A(k,q), are given a subscript according to the 
indices above with the number of cut propagators as a prefix. The double loop topology 
shown in figure 6(a) would be represented as, 

^s+t-l;ii---ir-i[ir,jt]ir+i---is-i[is,ji]j2---jt-i- 

In the case where the intersection vertex has no external legs attached, we represent it with 
a '*': 

As+t+3;ii---ir*ir+i---is*ji---jf i^-^) 

Following the same structure the butterfly topology, shown in fig. 6(b), would be repre- 
sented as, 

^s;ii---ir-i[ir,is,ji,j2]ir+i---is~i (D-5) 

We note that though the set of indices uniquely defines a topology, it does not account for 
possible symmetries between primitive amplitudes. The butterfly and higher multiplicity 
double loop topologies are of course beyond the scope of this paper. 
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